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ABSTRACT 
We report the discovery of the first eclipsing detached double white dwarf (WD) binary. In a pulsation search, 
the low-mass helium core WD NLTT 1 1748 was targeted for fast (w 1 minute) differential photometry with the 
Las Cumbres Observatory's Faulkes Telescope North. Rather than pulsations, we discovered w 180 s 3%-6% 
dips in the photometry. Subsequent radial velocity measurements of the primary white dwarf from the Keck 
telescope found variations with a semi-amplitude /^i = 27 1 ± 3 km s~' , and confirmed the dips as eclipses caused 
by an orbiting WD with a mass M2 = O.648-O.771M0 for Mi = O.1-O.2M0. We detect both the primary and 
secondary eclipses during the PqA = 5.64 hr orbit and measure the secondary's brightness to be 3.5% ±0.3% of 
the primary at SDSS-g'. Assuming that the secondary follows the mass-radius relation of a cold C/O WD and 
including the effects of microlensing in the binary, the primary eclipse yields a primary radius of Ri = 0.043- 
O.O397?0 for Ml = O.I-O.2M0, consistent with the theoretically expected values for a helium core WD with 
a thick, stably burning hydrogen envelope. Though nearby (at « 150 pc), the gravitational wave strain from 
NLTT 11748 is likely not adequate for direct detection by the Laser Interferometer Space Antenna. Future 
observational efforts will determine M\, yielding accurate WD mass-radius measurement of both components, 
as well as a clearer indication of the binary's fate once contact is reached. 
Subject headings: hinaiies: eclipsing — stars: individual (NLTT 11748) — white dwarfs 



1. INTRODUCTION 

Double white dwarfs (WDs) in tight enough (forb < 
day) binaries to reach contact in a Hubble time a re ex- 
pected on theoretical grounds (Nelem ans et al.l [2001), and 
are presumed to be the progenitors of highly variable ob- 
jects: R CrB stars, AM CVn binaries, and T ype la su- 
pernovae (llben & Tutukovlll984t IWebbink|[T98l . However, 
examples of thes e systems are rare, with only 10 known 
prior to our work (Nelemans et al. 2001; Badenes et al. '2009| 
MuUallv et al. 2009; Kilic et a l. 2009; Marsh etal. 201 



Kulkarni & van Kerkwiikll2010l) . Many of these binaries are 



of immediate interest for the Laser Interferometer Space An- 
tenna (LISA), provid ing 'verification' sources loud enough to 
be detected in space ( Nelemans '2009|) . 

Spectral measurements of the high proper motion object 
NLTT 11748 bv iKawka & Vennesi (2009) revealed it to be a 
low-mass (< O.2M0) helium core WD with logg = 6.2 ±0.15 
and Teff = 8540±50K and an H-rich (DA) atmosphere. The 
strong expectation that low-mass He WDs can only be form ed 
in common envelope events led iKawka & Vennesi (l2009h to 
note that this object was likely a binary. Indeed, they found a 
40km s~' velocity difference between two measured spectra. 
NLTT 1 1748 was originally targeted in a ZZ Ceti like pulsa- 
tion search based on theoretical calculations and the \ogg and 
Teff measurements (Steinfadt et al. 2010). In this failed pul- 
sation search, our rapid (« 1 minute) differential photometry 
with the Las Cumbres Observatory's Faulkes Telescope North 
(FTN) revealed w 180 s 3%-6% dips, which were confirmed 
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by Keck spectroscopy to be primary and secondary eclipses 
from a w 0.7 Mq faint C/O WD companion orbiting at 5.6 hr 

NLTT 11748 joins t he growing c l ass fe awka et all 
2006; Kilic et al.' '2007'; 'Mullallv et al.' '200^ Kilic et all 
2009; Marsh etal. 2010; Kulkarni & van Kerkwijk 201(J 
IBadenes et al.ll20 09) of WD binaries where the low-mass WD 
is a very low-mass (< O.2M0) helium core WD. These are 
made when the lower mass star overfills its Roche lobe at 
the base of th e red giant branch (Ivan Kerkwijk et all 119961: 
ICallanan et al.l [1998; Bassaetal. 2006a), triggering a com- 
mon envelope event that leads to an in-spiral of the now ex- 
posed low-mass He core to a shorter orbital period. Theoreti- 
cal work has found a dichotomy in the evolution of these low- 
mass He co re WDs ( Driebe et al]|1999HSerenelli et al.ll2002l; 
iPanei et an l2007) that dramatically impact their cooling rate. 
For masses >O.2M0 (metallicity dependent), the H envelope 
undergoes shell flashes that reduce its mass, eventually allow- 
ing the WD to cool rapidly. However, less massive (< O.2M0) 
WDs undergo stable H burning for gigayears, dramatically 
slowing their evolution and keeping them brighter for much 
longer than expected. 

Our observations are summarized in Section |2] demon- 
strating that the companion to NLTT 11748 is a cold (Teff < 
7400 K), old (1.5-3 Gyr) C/O WD. We show in Section[3]that 
eclipse modeling yields the first radius measurement of an ex- 
tremely low-mass WD. We close in Section |4]by discussing 
the value of additional measurements for constraining WD 
mass-radius relations. 

2. OBSERVATIONS 

We targeted NLTT 1 1748 as part of an observational search 
for ZZ Ceti like pulsations from low-mass He core WDs 
(ISteinfadt et al.ll2008bl I2OIOI) . Our photometric discovery led 
to our Keck spectroscopy, both described here. 

2.1. Faulkes Telescope North Photometry 
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TABLE 1 

Measured and Derived Quantities for NLTT 1 1748 
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Error Bar 
2009-12-22 
2010-01-08 
2010-02-03 



-0.2 0.0 



0.2 0.4 0.6 
Phase (cycles) 



0.8 



1.0 



1.2 



Fig. 1 . — Top panel: phased radial velocities from Keck/HIRES Ha mea- 
surements over two consecutive nights. The black line is the (e = 0) sinusoid 
with Ki = 271kms"' amplitude and the v,. = 133kms"' systemic velocity 
offset. Inset shows the phase combined Ha spectrum, in arbitrary flux units, 
with fitting function over plotted. Bottom panel: phased photometric light 
curve for all three nights of data. The eiTor bar represents 4 — 5 mmag uncer- 
tainty. 

Phot ometric o bservations were done with the 2 m robotic 
(FTN) jLewis e t al. 2010), part of the Las Cumbres Obser- 
vatory Global Telescope (LCOGT) networlQ, on Haleakala, 
Hawaii. We used the Merope camera with 2x2 pixel binning 
for a pixel scale of 0"28 pixel"' and a 4.75' x 4.75' field-of- 
view. Our observational setup, consisting of the SDSS-g' fil- 
ter and 45 s exposure time (with «22 s dead time), was aimed 
at detectin g pulsations, expected to be at the few minutes 
timescale ( Steinfadt et alj|2008bl 120101) . Preliminary reduc- 
tion, including bias and flat-field corrections using standard 
IRApQ routines, is done automatically at night's end. 

Our first observation on 2009 December 22 (night 1) was 
for 4 hr We immediately identified two dips in the light curve, 
«3 minutes in duration and a few percent in depth, indicative 
of eclipses. To verify the existence of eclipses we observed 
the system again on 2010 January 8 (night 2), for 3.5 hr, re- 
vealing two similar eclipses at the expected times based on 
the previous observations. An additional 30-minute observa- 
tion was obtained on 2010 February 3 (night 3). 

Photometric processing was carried out using IRAF aper- 
ture photometry and the VAPHOT package of routines 
dPeeg & Dovlell200l]) . Three comparison stars were selected 

'' http://lcogt.net 

' IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. http://iraf.noao.edu 



Quantity 


Value 


HIRES spectra 


Rad. vel. amp., Ki (kms"') 
Sys. radial velocity, iv (kms"') 
X^/dof 


271(3) 
133(2) 
4.5/6 


FTN photometiy 



Time of prim. eel. (BJD TDB) 


2,455,196.87828(7) 


Period (days) 




0.2350606(11) 




Ephemeris x /dof 




3.5/3 




Primary eclipse depth, di 




0.067(3) 




Secondary ecUpse depth, ^2 




0.034(2) 




Primary eel. duration, ri (s) 




180(6) 




Secondary eel. duration, ri (s) 




185(10) 




F2/F1 (in SDSS-g') 




0.035(3) 




Out of eclipse" X^^dof 




404.8/372 




Combined data (assuming 


Mass of Primary Mj = O.15M0) 


Limb darkening coefficient, hld 


0.0 


0.3 


0.5 


Mass of secondary, A/2 (Mq ) 


0.71(2) 


0.71(2) 


0.71(2) 


Inclination (deg) 


89.90(11) 


89.88(11) 


89.87(11) 


Radius of primary, S| (Rq) 


0.0393(9) 


0.0406(9) 


0.0415(9) 


X-/dof* 


285.1/227 


279.5/227 


276.5/227 


Distance*-' (pc) 




150(32) 




Sys. kin. (f/, y, B')''(kms-') 


(-151(9), -183(41), -. 


54(5)) 



Note. — Quantities in parentheses are la uncertainties on the last digit. 
" Excluding ±0.005 cycles around each eclipse. 
^ Result of fitting p hases ±0. 1 cycles a round each eclipse. 
'^ Scaled from Kawk a & VennesI #2009) using our Rj . 

'' Calculated from the pr oper mo tion ( Kaw ka & Vennesll2009l) corrected to 
the local standard of rest jHogg et al.l2005l) and our updated distance. 

to yield a light curve with the smallest variance. In 45 s ex- 
posures, each of these stars accumulated « 4xl0^e~ while 
NLTT 1 1748 accumulated w 7 x 10^ e". The differential fight 
curve for NLTT 1 1748 was constructed usin g the scheme de- 
tailed in Steinfadt et al. (2008b, inspired by iSokoloski et alJ 
I2OOII) . A de-trending second-order polynomial removed un- 
corrected long-term variations, such as color-airmass effects. 
For all data points outside of eclipses we find a x^/dof of 
404.8/372, indicating a reasonable estimate of our uncertain- 
ties. 

After considering the radial velocity measurements in 
Section 12.21 we found that nights 1 and 2 had primary 
and secondary eclipses, while night 3 shows a secondary 
eclipse. Using the technique of iKwee & van WoerdenI 
(1956) we fit the transit centers of each ecUpse, and 
combined all eclipse timings (primary (BJD TDB): 
2,455,188.886207(89), 2,455,205.81071(29); secondary 
(BJD TDB): 2,455,188.76882(16), 2,455,205.92783(19), 
2,455,231.784879(19)) in a linear least-squares analysis 
to find the orbital period of Porh = 5.641454 ±0.000026hr 
(barycentering used the JPL DE405 ephemeris). Fitting 
ephemerides using the primary and secondary eclipses sep- 
arately gives consistent results with our reported ephemeris 
(Table [T]l- The resulting phased light curve from all three 
nights of data is shown in Figure [1] 

2.2. Keck Spectroscopy 

We observed NLTT 11748 with Hig h Resolution EcheUe 
Spectrometer (HIRES, 'Vogt et al.l ll994 on the 10-m Keck I 
telescope 6 times on 2010 February 2 and twice on 2010 
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February 4. All observations were taken with the same grating 
tilt, integrating for 5 minutes, binning by 3 pixels in the spatial 
direction, and with a 0"86x 14" slit (giving /Jw 50, 000). The 
wavelength solution was based on a Th-Ar lamp, accurate 
to < 0.1 pixel (with 1. 4km s~' pixel"'), and the data covered 
3600- 8000 A. The stability of the solution was monitored to 
high precision by G. Marcy (private communication, 2010). 

We used the Ha line for our radial velocity measurements, 
although other lines were detected at lower significance. We 
initially identified the velocity from each spectrum by eye, 
then shifted each spectrum to zero velocity and combined 
them. We fit the combined Ha spectrum with the sum of 
a broad (« 520km s~' FWHM) Lorentzian and a narrow (« 
47kms~' FWHM) Gaussian with roughly equal depths. Us- 
ing this template, we then measured radial velocities and fit 
them to a sinusoid with the period and phase constrained by 
the eclipse timing. We then iterated this procedure, using the 
velocities from the radial-velocity fit to construct the template, 
fitting the template shape, fitting the individual velocities, and 
fitting the radial velocity curve. This converged quickly: after 
one iteration all changes were < 1 km s~' . 

We obtain a x^ for the radial velocity fit (assuming a circu- 
lar orbit and an accurate ephemeris from the eclipse photom- 
etry) of 4.5 for 6 dof. Our final radial-velocity amplitude was 
^1 = 271 ± 3kms~' yielding a mass function O.48(2)M0, or 
an M2 = 0.71{2)Mq secondary for a fiducial Mi = O.15M0 
primary. We computed the systemic radial velocity from 
the data, correcting to the solar system barycenter, finding 
Vr= 133 ±2kms~ . We do not correct for gravitational red- 
shift as uncertainties in Mi and Ri are too high and we expect 
a shift of only w 2 km s~' . We checked our velocities by fitting 
only the broad component of the Ha line or the H7 line, and 
both gave consistent results although with a factor of 6 lower 
precision. We limit the eccentricity of the orbit to < 0.06 at 

3(7. 

At the primary's orbital velocity, the width of the narrow 
Ha component may be broadened by orbital motion: the 
maximum derivative is InKi/Poib = 0.08kms~'s"^ Com- 
pared with the narrow Ha component, « 47km s~' FWHM, 
this effect is unlikely to affect our data (5-minute exposures) 
but might have slightly influenced the log^ measurement of 
iKawka & Velines (2009) . 

3. BINARY RA.RAMETER ANALYSIS AND RESULTS 

The deepest (primary) eclipse is the transit of the primary 
object (1) by the secondary (2), and the shallowest (sec- 
ondary) is the occultation of the secondary by the primary. 
There are six relevant binary parameters to be fitted: the pri- 
mary and secondary radii and masses (Ri, R2, Mi and M2); 
flux ratio (for our filter/CCD response) F2/F1; and inclina- 
tion /. We consider the orbital period, Porb, and the linear limb 
darkening coefficient, Mld. as fixed. We presun ie Mi ,n = 0.3 but 
will show ho w its uncertainty (mld ~ .2-0. 6, lLittlefair et al.l 
120061 120071 and iMaxtedetal] 120071) modifies our desired 
parameters. The eclipse light curve in Figure |2] shows the 
secondary eclipse as flat-bottomed while the primary eclipse 
shows the distinct signature of limb darkening. This indicates 
that the eclipses are likely total (confirmed by our inclination 
measurements) and R2 < Ri, F2/F1 ^ 1, and / w 90°. 

We initially fit the eclip ses with a parame terization of the 
linear limb darkening law dVan Hammdl 19931) for the primary 
ecHpse and a box function for the total secondary eclipse. 
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Fig. 2. — Phased light curves of the primary (top panel) and secondary 
(bottom panel) eclipses. The widths of each data point represent the exposure 
duration. The black line is the light curve model for the mld =0.3 parameters 
in Table[T] The crosses show the average of the model over the exposure du- 
ration of that associated data point. It is the differences between these points 
and those measurements that are considered when calculating x^ ■ Insets dis- 
play the la, 2a, and 3a contours of the x' minimization for the eclipse depth 
and width parameterized light curve fitting. The error bars are the [a error 
bars on the individual parameters. 



yielding the depths, di and d2, and durations, ri and T2. The 
1(7 results are reported in Table [T] with the la, la, and 3a 
contours plotted in the insets of Figure |2] This allows us to 
directly measure F2/F1 , as \-d2 = Fi/ {Fi+ F2). 

The sparsity of our data in eclipse requires additional as- 
sumptions about the properties of the secondary. The flat 
bottom secondary eclipse unambiguously gives a low flux ra- 
tio (at SDSS-g') of F2/F1 = 0.035(3) which implies that the 
value of the primary eclipse total depth must be dominated 
by the radii ratio. Our measured primary eclipse depth im- 
plies R2/R1 ~ 1/4, further corroborated by the unseen (and 
thereby fast) ingress and egress in the secondary eclipse. Our 
measurements of F2/F1 and R2/R1 constrain the tempera- 
ture of the secondary relative to the 8500 K primary. Taking 
SDSS-g' measurements as bolometric (L2/L1 =F2/Fi) gives 
T2 < 7400 K, which we confirrn by using the synthetic pho- 
tometry of 'Holberg & Be rgeronI (12006). constraining the flux 
ratio of model atmospher es in tegrated over the filter passband. 
The lKawka & VennesI (I2009h spectrum constrains the primary 
to be a low-mass, < 0.25 Mq, likely He WD with a radius 
« O.QARq. This implies an R2 « 0.01 Rq companion. There- 
fore, the secondary cann ot be a main sequence (MS) star as 
iKawka & VennesI (120091) constrain such an MS star to have 
M < Q.2Mq based on Two Micron All Sky Survey photom- 
etry, breaking our mass function and radius constraints. A 
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constant logg. 



brown dwarf would break our mass function and radius con- 
straints. A neutron star is far too small and its luminosity 
contribution insufficient. Therefore, a cold C/O WD is the 
only solution, meeting the radius, mass, and luminosity con- 
straints. 

We powerfully constrain the binary parameters by directly 
integrating the faces of the primary and secondary and model 
our light curve for any arbitrary inclination and phase by as- 
suming a linear limb darkening law. W e do this with a code 
originating from lSteinfadt et al I (l2008ah and Ma ndel & Agoll 
(|2()02|) . including the effects of microlensing in the binary 
(following 'Marsh 20011 and |^ol 2002). For the approximate 
parameters above, the Einstein radius is Re = ^J^GMjajcP- « 
0.003 7?0, giving RxJRe « 13 and RijR^ « 3, leading to a 
1 % magnification of the primary at conjunction. Our mea- 
surement QiK\ and FijFx reduces our six-parameter system 
to four-parameter system. Our surface temperature measure- 
ment yields an age of 1.5-3.0 Gyr (Chabri er et al.ll2000l) for 
the secondary allowing us to use M-R relations of cold WDs 
^Ithaus & Benvenuto 19981) for our current work. Though 
iKawka & Vennesi (i2009l) derived an My = 0.167 ± 0.005 Mg 
from their low-resolution spectrum, radial velocity broaden- 
ing and theoretical model uncertainties may render this deter- 
mination unreliable. Therefore, we calculate separate x^ grids 
of Ri and inclination for a range of Mi . Minimization of x^ 
on each grid constrains Ri and inclination directly and Ax^ 
techniques yield their uncertainties. We scale these uncertain- 
ties to consider the additional uncertainties of F2/F1 and Ki 
using a Monte Carlo method. 

The results of our direct eclipse fitting method are displayed 
in Figure |3] with the parameters and uncertainties reported in 
Table[T]forMi = O.ISMq. Without including lensing the fits 
would be very similar, but the inclination would be slightly 
lower (closer to 89.4°) and the radius of the primary would 
be slightly larger (0.0437/^q for mld = 0.3). The x^ min- 
ima for all Ml within each mld data set in Figure |3] are ap- 
proximately consistent and as such yield no constraint on Mi 



over this range. Our range of mld = 0.0, 0.3, and 0.5 results 
in a « 5% spread in Ri measurements, consistent with the 
log g = 6.2 ±0. 15 of lKawka& Vennesi (l2009b . The x^ value 
is reduced slightly as mld increases; however, we cannot con- 
strain Mld due to the sparsity of data at ingress/egress. This 
will be constrained in future higher cadence data sets. For 
Ml > Q.2Mq the minimum x^ values significantly increase 
and the quality of the eclipse fits degrade. Therefore, our 
combined radial velocity and light curve data sets provide an 
upper limit, assuming that the secondary is a C/O WD, of 
Mi<0.2Mq. 

4. CONCLUSIONS 

Our discovery of the first eclipsing detached double WD 
binary has substantially constrained low-mass He WD mod- 
els with rad ius measurements that agree with current mod- 
els (ISerenelli et al. 2002 : Panel et al . 2007). Additional mea- 
surements will reveal even more about the faint w 0.7 Mq 
C/O WD secondary. With our discovery, there are 1 1 con- 
firmed double WD binaries that will merge within a Hub- 
ble time. From spectroscopic masses alone, six of these 
contain low-mass (< 0.25Mq) He core primaries with large 
(> 0.03R(7)) radii indicative of a stable H burning shell 
dMuUallv et all 120091: iKiUc et al.ll2009t iBad enes etalJ 120091: 
iMarsh et all 120 lot iKulkarni & van KerkwijkHIoiOl). Other 
than its inclination, NLTT 11748 is not extraordinary in any 
way relative to these, many of which offer larger solid angles 
for eclipse observations than NLTT 1 1748. 

Our analysis of this system is far from complete, but our 
existing constraints, Table[T]and Figure[3] highlight the value 
of further studies. High cadence photometry, resolving the 
ingress/egress of each eclipse, would yield the component 
radii as functions of inclination and orbital semimajor axis. If 
future spectroscopic observations can reveal the secondary's 
line features and measure its radial velocity, this would allow 
precise measurements of the masses of both components. De- 
tection of eclipses in different filters, both bluer and into the 
infrared, would additionally constrain the temperature of the 
secondary further as well as its age. 

Model-independent mass and radius determinations 
for He and C/O WDs are rare. For He WDs, there 
exist several pulsar binary systems for which Shapiro 
time delay measurements have yiel ded accura te mass 
determinations: PSR B1855-H09 JKaspi et alj [1994), 
PSR J0437-4715 (van StrateneLalJ'2001), PSR J0751H-1807 
(Bassa et al . 2006b), PSR J1012H-5307 (Callanan et al. 199^ 
van Kerkwii k et al. 1996), PS R J1911-595 8A (Bassa eLMJ 
2006a|), and PSR J1909-3744 (iJacobv et al.ir2003 . 2005"t). Of 



these, only four have radii measurements, PSR J0437-4715, 
PSR J0751H-1807, PSR J191 1-5958A and PSR J1909-3744. 
However, all rely upon either distance measurements and 
evolutionary model-dependent luminosities or gravity mea- 
surements derived from atmos pheric models. Therefore , 
He WD evolutionar y models JAlthaus & Benvenutol [19981: 
ISerenelU et akl 120021: iPanei et al.l l2007h remain largely un- 
constrained. For C/O WDs, mass and radius measurements 
have been more frequent (SchmidJ Il996t [Provencal et al.l 
119981 120021: Eisewell et al. 2009) and all rely on distance 
measurements via parallax or cluster membership, surface 
brightness measurements and detailed knowledge of each 
spectrum. To reach the precisions required to constrain 
theoretical models, these measurements must rely upon 
atmospheric model spectra. Only recently Parsons et al.l 
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( 120 lOl) have derived precise measurements of mass and radius 
in a model-independent way in the eclipsing binary system 
NN Ser NLTT 11748 offers an additional C/O WD system 
for model-independent mass and radius measurements as 
well as the first such system for He WDs. 

Gravitational wave emission will bring this system into con- 
tact in 13.8-6.3 Gyr (for Mi = 0.1-0.2Mo). With a mass 
ratio of 0.15/0.71 = 0.2, stable mass transfer a t contact is 
likely to occur and create an AM C Vn binary ( Nel emans et al.l 
1200 U iMarsh et al.l 12004 . The known orbital period, in- 
clination, and distance (d « 150pc) yield th e gravitational 
wave strain at Earth, h = (3.6-7.9) x 10"^^ (iTimpano et al.l 
l2006t iRoelofs et al.l l2007h at a frequency z/ « lO'^^Hz for 
Ml = Q.1-Q.2Mq. Though not 'loude r' than the verification 
sources tabulated bv lNeleman"sl (l2009h . NLTT 11748's accu- 
rate ephemeris allows for a coherent folding of the LISA data 
over the mission duration. 

We thank the referee for helpful comments and suggestions. 
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